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Abstract 

In several recent papers it was claimed that SN 1987A data supports the existence of 4.0 eV 
and 21.4 eV active neutrino mass eigenstates, and it was suggested that such large active neutrino 
masses could be made consistent with existing constraints including neutrino oscillation data and 
upper limits on the neutrino flavor state masses. The requirement was that there exist a pair 
of sterile neutrino mass states nearly degenerate with the active ones, plus a third active-sterile 
doublet that is tachyonic (m 2 < 0). Here, independent evidence is presented for the existence of 
sterile neutrinos with the previously claimed masses based on fits to the dark matter distributions 
in the Milky Way galaxy and four clusters of galaxies. The fits are in excellent agreement with 
observations within the uncertainties of the masses. In addition, sterile neutrinos having the sug- 
gested masses address the "cusp" problem and the missing satellites problem, as well as that of the 
"top down" scenario of structure formation - previously a chief drawback of HDM particles. Nev- 
ertheless, due to the highly controversial nature of the claim, and the need for two free parameters 
in the dark matter fits, additional confirming evidence will be required before it can be considered 
proven. 
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I. INTRODUCTION 



In several recent papers evidence was presented based on an analysis of SN 1987A for the 
existence of two active neutrino mass eigenstates having masses 4.0±0.5e\^ and 21.4±1.2eV, 
hereafter referred to as the "claim" [llS]. The present paper reviews that claim and presents 
new independent evidence supporting it. It is shown that the problem of dark matter 
within galaxies and galaxy clusters can be resolved using a pair of sterile neutrinos having 
the very same masses as claimed previously. Even though the claim was based on active 
neutrinos from SN 1987A, their existence would require sterile neutrinos having nearly the 
same masses, so as to be compatible with various constraints. We first discuss the model 
for the sterile neutrino fits to the dark matter distribution in the Milky Way and galaxy 
clusters, then briefly review the prior work making the claim of 4.0±0.5el / and 21.4±1.2eV 
neutrino masses, and finally discuss some ways the claim could be corroborated further. 



II. STERILE NEUTRINOS AS HOT DARK MATTER IN GALAXIES AND 
CLUSTERS 



Traditionally, the dark matter problem has been mainly addressed by the assumption it is 
cold. However, recent literature indicates that cold dark matter cannot explain the existence 
of flat cores at the center of clusters 3^(| and the relative absence of satellite galaxies^]. 
The existence of warm or hot dark matter better reconciles theory and the observations by 
solving these problems. In fact, simulations show that a flat core can be produced and the 
number of dwarf galaxies surrounding the Milky Way would be smaller if the dark matter 
particles are warmer . Furthermore, since the sterile neutrino mass used is larger in 
the galaxies than in clusters, the formation time of the galactic halo will be shorter. Hence, 
in this scenario, the galaxies would completely form before the clusters, which matches the 
bottom-up formation scenario in the structure formation. If they exist most of the sterile 
neutrinos can be produced through active-sterile neutrino oscillation by resonant or non- 
resonant mechanisms during the early epoch of the Big Bang 10|, yj|. If sterile neutrinos 



are the dark matter particles, their halo profile should match the dark matter profile in both 
clusters and galaxies. 
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A. Derivation of sterile neutrino dark matter mass profile 



In this section, we derive the sterile neutrino halo profile and compare it to the observed 
mass profile in clusters and the Milky Way. 

The pressure of the sterile neutrinos is given by 
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where m s , g s , T s , e = p 2 /2m s , /i and p are the rest mass, degrees of freedom, temperature, 
kinetic energy, chemical potential and the momentum of the sterile neutrinos respectively 
If m s ~ 1 eV for a 1O 15 M cluster mass and m s ~ 20 eV for a 1O 11 M galactic mass, 
the number densities are about 10 5 cm" 3 and 10 6 cm -3 respectively. For T s ~ 1 K, this 
number density is nearly the same order of magnitude of the quantum number density 
hq = (2m s kT s /h 2 ) 3 ^ 2 ~ 10 6 — 10 7 cm" 3 . Therefore, most sterile neutrinos are in the classical 
regime or they are only slightly degenerate (e* 5 ^ 3> 1). Then, Eq. (1) can be written as 
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Let u = ^Je/kT s = pj \ / 2m s kT s and v = ^2e/kT s = pj \/m s kT s , Eq. (2) becomes 
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The integrals in Eq. (3) are both equal to 3yfn/8. Therefore we get 
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Additionally, the number density of the sterile neutrinos is given by 
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By using the same approximation method, we get 
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In the classical limit e^ kTs — > 0, Eq. (7) reduces to the ideal gas law P = nkT s . However, if 
we expand e^ kTs up to first order to allow a slight degeneracy, we have 



n 
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So that from Eq. (6), we then have 
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Substituting Eq. (9) into Eq. (8), and writing n as n(r) to reflect its dependence on radial 
distance, we obtain for the final equation of state for the slightly degenerate neutrinos 



P(r) = n(r)kT s + ^^( ms kT s )- s / 2 (kT s )-^ 2 . 

By using the hydrostatic equilibrium condition, since p{r) = m s n(r) we find 

dP(r) Gm s M(r)n(r) 
dr r 2 

and finally, the overall mass profile M(r) can be found as the sum of the dark matter profile 
and the observable baryonic mass profile M B {r) 
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M(r) = M D (r) + M B {r) 
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The three coupled eqs. (10-12) cannot be solved in closed form, but a numerical solution is 
possible if one assumes values for the parameters m s , T and po an d then iterates in radial 
steps outwards from r = 0. 

This leaves only the issue of the function Mr (r) . In clusters, the baryonic mass profile is 



commonly expressed in the /3-model profile 
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where pso, r c and /3 are the parameters that can be fitted by the surface brightness profile 
of a cluster. In the Milky Way, the baryonic mass is mainly composed by the bulge and 



disk. These components can be modelled by the gravitational potential model 
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where M cn , Md n , r cn , c n and b are fitted parameters in the model [13f| . Since the sterile 
neutrino temperature is nearly homogeneous, just like the cosmic microwave background, 



we can assume that the sterile neutrino is isothermal in an entire cluster and in the Milky 
Way. 

Note that the larger the sterile neutrino mass, the smaller will be the halo size. Therefore 
it is assumed that the dark matter profile in the Milky Way is mainly characterized by the 
21. 4 ±1.2 eV sterile neutrinos and that of galaxy clusters by the 4.0 ±0.5 eV sterile neutrinos. 
In Fig. 1 and 2, we fit the derived total mass profiles (from Eqs. 10-14) together with the 



data from four large c 
of the Milky Way [l5 



usters (A478, A1413, A2204 and PKS0745) |l4j and the rotation curve 
respectively. The central density of the sterile neutrino halo p s0 



and the sterile neutrino temperature T s are the only free parameters adjusted in the fits. 
The neutrino temperatures used in Figs. 1 and 2 are T s = 0.1 — 0.5 K (for m s = 4.0 eV) 
and T s = 0.07 K (for m s = 21.4 eV) respectively. The parameters used in the baryonic mass 



profile from each cluster can be obtained from observations 17H23I] . As can be seen, the fits 
agree very well with the observed data, suggesting that 4.0 eV sterile neutrinos may be in 
fact the major component of dark matter in galaxy clusters within 1.5 Mpc, and that 21.4 
eV sterile neutrinos may be its main component in the Milky Way within 25 kpc. 



III. PRIOR WORK MAKING THE CLAIM 



The values of the neutrino masses have been a matter of considerable interest both 
theoretically and experimentally ever since they were shown to be nonzero by virtue of the 
existence of neutrino oscillations. Empirically, for the flavor state masses only upper limits 
exist from various decay processes, such as tritum beta decay, wh.ch as of 20i3 yields an 
upper limit on the electron neutrino mass m < 2.0eV [24j. An even lower limit exists on the 
sum of the 3 flavor state masses from cosmology which can be as low as 0.24eV, depending on 
how the analysis was done [24| . It is of course only the mass differences between pairs of the 
mass eigenstates, Am 2 , that are revealed in neutrino oscillation experiments not the masses 
themselves, with the relation between the flavor and mass eigenstates and their respective 
masses being through the mixing matrix [Z^: 

(a) v Fi = ZUijVMj (b) m 2 Fi = T\Uij\ 2 rr? Mj (15) 

The ability of an experiment to measure the mass itself for a mass eigenstate is severely 
limited for eV-scale masses, and perhaps the only feasible sources for doing so are the 



5 



supernovae for which the long baseline and briefness of the initial pulse offer a hope of 
having the arrival times of distinct mass states separate out en route to Earth. The recent 



pair of papers making the claim of 4.0 eV and 21.4 eV masses was based on an ana! 
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ysis 



26). 



of SN 1987A 

HQ 

- and relies on an idea that was first investigated much earlier 
The recent papers, however, extended the earlier analyses by reducing the uncertainties in 
the claimed masses, and showing a way such large neutrino masses could be compatible 
with what is empirically known about neutrino mass limits . The claim rests on the 
assumption that most of the 24 neutrinos observed from SN 1987A were emitted during 
a time much less than the 15 second burst itself, possibly in as short a time as 1 second. 
If one simplistically assumes simultaneous neutrino emissions then the individual neutrino 
arrival times at the detectors t can be expressed in terms of the light travel time from the 
supernova to and that neutrino's mass m and energy E. Based on relativistic kinematics 
1 - m 2 c 2 /E 2 = v 2 /c 2 = t 2 /(t + tf « 1 - 2t/t , so that 

E 2 = ^ 
Note that by Eq. (15), t = is the time that infinite energy (v = c) neutrinos would 

have arrived at the detector - which given the neutrino time distribution in the detectors, 

is probably within ±0.5s the same as setting t — for the first arriving neutrino in each of 



the three unsynchronized detectors jl]. Note that Eq. (15) implies that on a plot of 1/E 2 
versus t all neutrinos having a fixed mass m lie on a straight line of slope M = 2/(m 2 io)- 
Thus, if neutrinos having two distinct masses were present among the 24 they would all lie 
close to one of two straight lines having slopes Mi = 2/ (mfto) and M 2 = 2/(m|to) that pass 
through the origin. Remarkably, as shown in, this pattern is exactly what the data show 



ay 



26]- 



(within the uncertainties in 1/E ) for every one of the 24 observed neutrinos 

There are many reasons for skepticism about the claim. One reason is that it is unlikely 
that all the 24 SN 1987A neutrinos were in fact emitted during a time as short as Is, and some 
core collapse modelers calculate that much of the neutrino emission occurs over an extended 
;ime. 27|] Nevertheless, other modelers do show the emissions to be very concentrated in time 
28]. For example, (28] shows 73% emitted in the first second and 82% in the first 2 s (values 
based on a numerical integration of their Fig. 1). Finally, we note that the goodness of the 
two straight line fit (from which the claim was inferred) is not significantly impaired if there 
is a limited spread in emission times. Thus, in a Monte Carlo simulation in which 24 events 



are generated, assuming an emission time distribution and energy spectrum following those 



of 



28], and assuming the masses of individual neutrinos are as claimed, the agreement of 



the fit to two straight lines is often just as good as for the actual data. 

Achieving consistency between such large claimed masses and existing constraints on the 
neutrino masses requires an unconventional model of 3 active and 3 sterile neutrino mass 
states. In this 3 + 3 model there are three nearly degenerate active-sterile doublets, as 
depicted in Fig. 1 of |2j. The masses of the first two active-sterile doublets are mi = 4.0eU 
and vri2 = 21AeV and the third active-sterile doublet mass is tachyonic, and has a conjec- 
tured mass: m| = — 0.2keV 2 . The splitting of the members of the two tardyonic (m 2 > 0) 
doublets are the values from neutrino oscillation experiments for Am^, and Am 2 tm respec- 
tively, while the splitting for the tachyonic doublet is taken as leV 2 - the Am 2 suggested 
by short baseline reactor neutrino experiments [29]. These choices make the fractional mass 
splitting of each doublet, ^p-, identical. It has been discussed elsewhere how this 3 + 3 
model could be consistent with various observations about neutrinos, including upper limits 
on the flavor state masses, neutrino oscillation data, and the upper limits on any tachyonic 
neutrino masses that might be the source of dark energy 3). Clearly, given its highly contro- 
versial nature the claim needs to be supported by other observations before it will be taken 
seriously. In fact, it is more impressive if such a confirmation comes from an entirely differ- 
ent area from supernova neutrinos. Therefore it is particularly interesting that the analysis 
reported in the present paper on the implications of sterile neutrinos for the dark matter in 
galaxies and galaxy clusters does in fact provide that independent sort of evidence. Unlike 
dark energy which depends on the average energy density throughout all space, dark matter 
tends to be concentrated in bound systems such as galaxies and clusters, and it is sensitive 
to the masses of the sterile neutrino mass eigenstates, and not the flavor eigenstates. 



IV. DISCUSSION 



In summary, here we have found excellent agreement with observations when using the 
claimed masses of two sterile neutrinos, 4.0±0.5eU and 21.4±1.2eU that had been inferred 
previously from an analysis of SN 1987A. Although such good fits (all within the cited 
la errors on the masses) provides confirmation for the claimed masses existence, the highly 
controversial nature of the claim and the use of two free parameters to achieve the fits means 
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that additional observations will be needed to prove the claim. The most direct confirmation 
might involve a neutrino oscillation experiment sensitive to Am 2 = 21. 4 2 — 4.0 2 = 442eV 2 . 
For example, at a 30 MeV neutrino energy we find an oscillation wavelength of about 0.167 
m, which could be readily observed, particularly using a high intensity source with a strong 
30%) monochromatic component, such as the Oak Ridge Spallation Neutron Source (SNS) 



30] . If we are fortunate enough to have another supernova in our galaxy any time soon, 



the good sensitivity of Super-Kamiokande and other detectors also could give such a clear 
confirmation or refutation of the claim. In fact, Super-K is even sensitive enough to search for 
a supernova burst from Andromeda (2 events expected). Surprisingly, even the observation 
of only 2-3 events associated with Andromeda could determine whether the claim of two 
mass eigenstates is correct, in light of the pulse spreading with source distance required by 
Eq. (15). However, it is important that such a search would need to be conducted differently 
from the one actually employed in a 2007 Super-K supernova search by Ikeda et al.[3l|, as 
described in a recent paper [321 ] . 
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FIG. 1: Observed and fitted total mass profiles for clusters A478, A1413, A2204 and PKS0745-191. 
Given the la errors on mi = 4.0 ± 0.5eV, the sterile neutrino masses used in the fits are 3.5 eV 
(dotted line), 4 eV (solid line) and 4.5 eV (dotted line) in each cluster. The central sterile neutrino 
densities and the sterile neutrino temperatures for the four clusters are: (p S Q,T s ) are (1.4 x 10~ 25 
g cm- 3 , 0.5 K) (A478), (1.5 x 10~ 25 g cm" 3 , 0.1 K) (A1413), (1.5 x 10~ 25 g cm" 3 , 0.5 K) (A2204) 
and (2.2 x 10~ 25 g cm" 3 , 0.4 K) (PKS0745-191) respectively. 
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FIG. 2: Observed and fitted rotation curve of the Milky Way. The solid dark curve is the observed 
data obtained by Sofue et al.(1999) and the error bars are the latest data obtained by Sofue(2011). 
The components of the sterile neutrino halo, disk and bulge that contribute to the total rotation 
curve are shown. The dashed line and the dotted lines are the rotation curves fitted by using 
m s = 21.4 eV (central density = 7.3 x 10~ 24 g cm" 3 ) and m s = 20.2(22.6) eV, (central density 
= 7.3 x 1CT 24 g cm" 3 for m s = 21 A eV and central density =6 x 1CT 24 g cm" 3 and 8.5 x 10" 24 g 
cm" 3 for m s = 20.2 eV and 22.6 eV respectively)The sterile neutrino temperature is 0.07 K. 
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